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We present the results of neutron-scattering studies on various aspects of crystalline and magnetic 
structure in single crystals of Lai.6- :E Ndo.4Sr a ;Cu04 with x = 0.12 and 0.15. In particular, we 
have reexamined the degree of stripe order in an x = 0.12 sample. Measurements of the width for 
an elastic magnetic peak show that it saturates at a finite value below 30 K, corresponding to a 
spin-spin correlation length of 200 A. A model calculation indicates that the differing widths of 
magnetic and (previously reported) charge-order peaks, together with the lack of commensurability, 
can be consistently explained by disorder in the stripe spacing. Above 30 K (i.e., above the point 
at which a recent muon spin-rotation study has found a loss of static magnetic order), the width 
of the nominally elastic signal begins to increase. Interpreting the signal as critical scattering from 
slowly fluctuating spins, the temperature dependence of the width is consistent with renormalized 
classical behavior of a 2-dimensional anisotropic Heisenberg antiferromagnet. Inelastic scattering 
measurements show that incommensurate spin excitations survive at and above 50 K, where the 
elastic signal is neglible. Given that the stripe order is believed to be pinned by the low-temperature 
tetragonal (LTT) crystal structure, we have also investigated the transition near 70 K from the 
low-temperature orthorhombic (LTO) structure. We show that our x — 0.12 crystal passes through 
an intervening less-orthorhombic phase, before reaching the LTT at ~ 40 K, whereas the x = 0.15 
crystal goes directly from LTO to LTT, with coexistence of the two phases over a range of ~ 7 K. 
Sharp Bragg peaks in the LTT phase of the x = 0.15 crystal indicate a domain size of > 1000 A, with 
no obvious evidence for LTO domains; hence, the coexistence of stripe order and superconductivity 
in this sample cannot be explained by a mixture of crystalline phases. Finally, we present scattering 
evidence for small LTT-like domains in the LTO phase of the x — 0.15 sample. A correlation 
between the volume fraction of such domains and deviations of in-plane resistivity from linear T 
dependence suggest that charge stripes interact with these domains within the LTO matrix. 

75.50.Ee, 75.30.Fv, 71.45.Lr, 71.27+a 



I. INTRODUCTION 

Important insights into the spatial correlations of 
charge and spin in cuprate superconductors have been 
provided by studies of the x = | anomaly originally 
discovered!!! in La2_ a; Ba 2 ;Cu04. Early on it was demon- 
strated by X-ray diffraction experiments that the anoma- 
lous suppression of superconductivity near x = g is asso- 
ciated with a modification of the crystal structure from 
the usual low-temperature orthorhombic (LTO) .-phase 
to the low-temperature tetragonal (LTT) variant nu Af- 
ter it was found that the same structural modification 
can be induced in. La2- a; Sr a; Cu04 by partial substitu- 
tion of Nd for Lap intensive efforts were focused on the 
latter system. A systematic study of the phase dia- 
gram revealed that not only the symmetry of the lat- 
tice modulation but also its amplitude is correlated with 
suppression of superconductivity and the appearance 
of local, stetic magnetism.El The discovery by neutron 
diffractionBQ of elastic superlattice peaks corresponding 



to two-dimensional charge and spin order in a sample of 
Lai.4gNdo.4Sro.i2Cu04 suggested a likely explanation for 
the anomaly: the dopant-induced charge carriers, whose 
density is spatially modulated in a periodic fashion sim- 
ilar to an array of stripes, can be pinned by the lattice 
modulation in the LTT phase, with the pinning being 
strongest when the periodicity is commensurate, near 
x = |. With the charge stripes localized, the Cu mo- 
ments in intervening regions can order antiferromagneti- 
cally, but with neighboring domains having an antiphase 
relationship induced by the charge. The suppression of 
superconductivity by the stripe pinning is compatible 
with at least one theoretical model for the cupratesja in 
which the development of superconducting phase coher- 
ence is limited by the ability of hole pairs to tunnel be- 
tween stripes. The Josephson coupling between stripes 
should be depressed by pinning. 

There is a clear connection between the static corre- 
lations found in h&i g_j.Ndo.4Sr3.CuO4 and the dynamic 
ones observed in good superconductors. For example, the 
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magnetic fluctuations in superconducting La 2 - a: Sr :l: Cu04 
are characterized .-by a doping-dependent incommensu- 
rate wave vectorja and, for a given Sr concentration, 
the wave vectors found i&,samples with and without Nd 
are essentially identicaLllj'til Also, elastic incommensu- 
rate magnetic peaks have now been observed in other 
superconducting systems, such as La2- a Sr :r Cui_j ; Zn a 04 



(Ref. ||, H, U) and La 2 Cu0 4+ «5 (Ref. |l5j). Furthermore, 
it has recently been shown that the magnetic fluctua- 
tions in underdoped,r^Ba2Cu306+x have an incommen- 
surate component JIS'EJ and that the spatial orientation 
of the modulation wave vector is identical with that in 
La2-sSr s Cu04. Thus, it may be possible to gain a better 
understanding of the dynamic correlations in the good su- 
perconductors by studying the static correlations in the 
Nd-doped system with depressed superconductivity. 

A number of intriguing features have already been ob- 
served in the stripe-ordered phases. Contrary to ob- 
servations on Fermi-surface-driven charge-density-wave 
(CDW) ordered systems,t3 infrared reflectivity studies 
indicate that there is no significant rgp-P in the optical 
conductivity within the CuC>2 planesJla Also, there is in- 
creasing evidence that buJJt-suaer conductivity and static 
stripe order can coexist OEj£il£3 Given these surprising 
results, it is important to investigate further the nature 
of the static stripe correlations. What is the nature of 
the order, both of the spins and of the underlying lattice 
modulation? 

We have previously shown that the spin and charge 
order in Lai^-zMAj^Sr^CuG^ with x = 0.12 is quasi- 
two-dimensionalJ3'E3 Here we show that the ordering is 
somewhat glassy. Using high Q resolution measurements 
of an elastic magnetic peak, we find that the peak width 
k saturates below ~ 30 K, at a value corresponding to a 
spin-spin correlation length (= l/«) of ~ 200 A. A model 
calculation, presented in the discussion section, demon- 
strates,— that the differing spin and charge-order peak 
widthsES, together with the deviation from commensu- 
rability, are consistent with disorder in the charge-stripe 
spacing. The cause of the disorder is undetermined, but 
might be due to interactions between the holes and the 
Sr 2+ dopant potentials. ._. 

A recent muon-spin-rotation (/iSR) studyE3 indicates 
that, for x — 0.12, the static magnetic correlations disap- 
pear at approximately 30 K. In contrast, the nominally 
elastic magnetic peak detected with neutrons does not 
disappear until ~ 50 K; however, we find that the peak 
width begins to grow for T > 30 K. The /iSR result im- 
plies that the higher temperature neutron signal involves 
an integration over low frequency spin fluctuations. Un- 
der the assumption that the finite energy resolution in- 
tegrates over the dominant critical fluctuations, k cor- 
responds to the inverse of the instantaneous spin-spin 
correlation length. The temperature dependence of k is 
consistent with the renormalized-classical regime of an 
anisotropic Heisenberg antiferromagnet.Ou A fit to the 
data yields an estimate of the anisotropy of the effective 
magnetic exchange; however, in order to obtain consis- 



tency with the narrow Q-width observed for spin fluctu- 
ations at Tilo = 3 meV, it is necessary to consider quan- 
tum effects such as those predicted when the magnetic 
domains correspond to even-legged spin ladders. lj'lZI 

We have also made a detailed study of the LTO to 
LTT transition in crystals with x — 0.12 and 0.15. One 
significant result is that, for the x = 0.15 sample, the 
LTT phase exhibits sharp Bragg peaks, with no obvious 
evidence for LTO domains. It followsJJxat the stripe or- 
der and superconductivity observedO'E3t3 in this sample 
must come from the same crystallographic phase. In con- 
trast, we do observed diffuse scattering consistent with 
LTT-like domains within the LTO phase for T > 80 K. 
The volume fraction of such domains is correlated with 
a deviation of the in-plane resistivity from linear T de- 
pendence. This correlation suggests that charge stripes 
may interact with the LTT-like domains, thus affecting 
the resistivity. 

The rest of the paper is organized as follows. After 
a brief description of experimental details in the next 
section, characterizations of the LTO-to-LTT structural 
transition are presented in section III. In section IV, the 
high Q resolution study of elastic magnetic scattering for 
x = 0.12 is reported. Some inelastic magnetic scattering 
measurements for the x=0.12 and 0.15 compositions are 
presented in section V. In section VI, elastic diffuse scat- 
tering measurements near a forbidden superlattice posi- 
tion in the LTO phase of the x = 0.15 sample are an- 
alyzed. Further discussion of the results, including the 
stripe-disorder model calculation, is presented in section 
VII. 



II. EXPERIMENTAL DETAILS AND NOTATION 

The crystals of Lai.g_ :E Ndo.4Sr a ;Cu04 used in this work 
were grown at the University of Tokyo by the traveling- 
solvent floating-zone method. All are cylindrical rods 
4 mm in diameter, with varying lengths. The elastic 
scattering studies of the x = 0.12 composition were per- 
formed on crystal U2, which is approximately 10 mm in 
length. The corresponding inelastic measurements were 
done on crystal U3, which was originally 40 mm in length, 
but which eventually fell into pieces due to hydration of 
small, poorly reacted inclusions. Pieces of the latter crys- 
tal were used in the X-ray diffraction study of charge- 
order scattering reported in Ref. For the x = 0.15 
composition, crystal U6 (35 mm in length) was used. An 
earlier crystal of the same composition (U4) was used for 
the measurements reported in Ref. [□]. That crystal suf- 
fered the same fate as U3, and a piece of it was later used 
for the magnetization study described in Ref. [|(]. 

The neutron-scattering studies were performed on 
triple-axis spectrometers H4M and H7 at the High Flux 
Beam Reactor, located at Brookhaven National Labo- 
ratory. Pyrolytic graphite (PG) monochromators and 
analyzers were used, and the (002) reflection was em- 
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ployed in all cases, except for certain elastic, high (ir- 
resolution measurements of strong Bragg peaks, where 
the analyzer was set to the (004) reflection. Essentially 
all measurements were done with an incident neutron 
energy of 14.7 meV, with one or more PG filters in the 
incident beam to minimize the neutron flux at higher- 
harmonic wavelengths. Each crystal was cooled in a Dis- 
plex closed-cycle He refrigerator, with the temperature 
monitored by a Si diode. 

As mentioned in the introduction, the three crystal 
structures of relevance to the present work are the LTO, 
LTLO, and LTT, which correspond to thc—Bmafr, Pccn, 
and P4 2 /ncm space groups, respectively.Qca All of these 
structures have a unit cell with lattice vectors a and b ro- 
tated by 45° with respect to the Cu-0 bonds in the CuC>2 
planes. When we specify a reflection in reciprocal lattice 
units (2ir/a , 2ir/b , 2tt/c) with respect to this cell, we 
will denote it with a subscript o for "orthorhombic" . (In 
the LTT phase at 10 K, a = 5.34 Aand c = 13.1 A.) To 
describe the magnetic scattering, it is more convenient 
to work with respect to a smaller cell with lattice vectors 
parallel to the in-plane Cu-0 bonds. In this case, the re- 
ciprocal lattice units are {2\/2ir/a , 2y/2n/b , 2n/c), and 
no subscript will be appended to the reciprocal-space co- 
ordinates. 



III. LTO-TO-LTT TRANSITION 

In order to interpret the nature of the magnetic order, 
it is first necessary to characterize the underlying crys- 
talline order. As a starting point, Fig. [j] shows the dif- 
ference between the in-plane lattice parameters, b — a, as 
a function of temperature for two samples. The x = 0.15 
crystal transforms directly from LTO to LTT, with a co- 
existence of significant fractions of the two phases be- 
tween 71 and 77 K (see Fig. ^| below). In contrast, the 
x = 0.12 sample shows a large but incomplete decrease 
in the orthorhombic splitting at approximately 69 K. 

The measurements on the x = 0.12 sample require a 
bit of explanation. In the LTO phase, there are four pos- 
sible orientations of twin domains. With the large mosaic 
crystals commonly studied, one typically finds the simul- 
taneous presence of all four domain orientations, with 
comparable populations of the four domains. As a result, 
a scan along the [100] direction through a Bragg peak at 
(H00) o will also detect a twin-related peak corresponding 
to (0-ff0) o . Such is the case for the x = 0.15 sample, and 
it allows one to measure the orthorhombic splitting with 
a single scan. With the x = 0.12 crystal (U2), on the 
other hand, a single domain dominates (and there are a 
total of 3, rather than 4, domains, with unusual relative 
orientations). As a result, it is necessary to separately 
scan the (H0Q) o and (0H0) o reflections. An advantage 
of this situation is that, because the peaks do not overlap, 
we can achieve a better resolution of small changes in or- 
thorhombicity. On the other hand, our absolute measure 
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FIG. 1. Temperature dependence of the difference in 
in-plane lattice parameters, b — a (based on LTO unit cell), 
for the Lai.6- :E Ndo.4Sr :E Cu04 crystals with x = 0.12 (circles) 
and x = 0.15 (squares). 



of orthorhombicity is limited by the precision with which 
the crystal is aligned. 

The splitting measured for the x = 0.12 crystal below 
the structural transition temperature is shown in Fig. |^. 
The error bar indicates twice the maximum variation in 
splitting between cooling and warming cycles (and cor- 
responds to 4% of the peak width). The absolute uncer- 
tainty is certainly much greater than this, and we cannot 
establish from these measurements whether the crystal is 
tetragonal or orthorhombic at low temperature. It is cer- 
tain that any remaining orthorhombicity below ~ 40 K is 
extremely small, and, based on powder diff taction stud- 
ies of samples with the same composition,ErErnia we ex- 
pect that the crystal is actually tetragonal. The signifi- 
cant point is that the splitting grows monotonically be- 
tween 50 K and the transition at 69 K. This behavior is 
consistent with the occurrence of an intervening LTLO 
phase, as first observed by Crawford et ala. A corre- 
sponding feature is the increase in integrated intensity 
for the (200) o and (020) o Bragg peaks (see Fig. ||), which 
was seen previously in a single-crystal study by Shamoto 
et ciIlQ The most likely explanation for the intensity vari- 
ation is a reduction in extinction due to strain associated 
with the growing orthorhombicity. Assuming this to be 
the correct interpretation, the intensity provides a more 
precise measure of the structural changes, and it suggests 
that the crystal is essentially tetragonal, and no longer 
changing, for T < 40 K. 

Turning to the x = 0.15 crystal, Fig. || shows high- 
resolution elastic scans through the (040) o /(400) o Bragg 
peaks. At 80 K, the LTO peaks are sharp and well sepa- 
rated; however, close inspection reveals a weak and broad 
contribution centered midway between the strong peaks. 
A similar feature has recently been observed in an_X-ray 
powder diffraction study by Moodenbaugh et al.pB and 
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FIG. 2. Comparison between the orthorhombic splitting, 
b — a, (circles) and the average intensity of the (200) and 
(020) Bragg peaks (squares), normalized at low temperature. 



interpreted as evidence for small LTT-like domains. We 
will return to that issue in section VI. For now, we fo- 
cus on the sharp LTT peak that begins to appear below 
78 K. As indicated in (b) and (c), there is coexistence of 
the LTO and LTT phases down to approximately 71 K. 
(Moodenbaugh et a/.cj observed a broader coexistence re- 
gion of approximately 30 K in their powder study.) While 
a single, strong LTT peak dominates at 70 K, there is also 
a very weak component with a width roughly 3 times 
greater that lies under it. 

The point to be emphasized here is the narrow width 
of the dominant LTT peak at 70 K. It is only marginally 
broader than the LTO peaks at 80 K, and all are close 
to the resolution width. Taking the resolution into ac- 
count, we estimate a typical domain size of - 1000 A or 
greater. The LTT peak observed in this crystal is con- 
siderably-,sliarper than those typically detected in powder 
sampleSjCraH where the LTT peak width is comparable 
to the LTO splitting. The broad peaks observed in pow- 
ders raised concerns that the superaanduGitivity which ap- 
pears to survive in the LTT phasec3'EllE3 might actually 
be associated with small LTO domains. In light of the 
present results, such a possibility seems quite unlikely. 



IV. ELASTIC MAGNETIC SCATTERING 

As discussed elsewhere IIQ magnetic superlattice peaks 
are observed within the (h, k, 0) zone of reciprocal space 
below approximately 50 K for x = 0.12, and 45 K for 
x = 0.15. Below ~ 3 K, where the Nd moments start to 
order via coupling to the ordered Cu moments, thus en- 
hancing the superlattice intensity with their much larger 
moments, it was possible to investigate the correlation 
lengths. Within the planes, the correlation length was 
found to be substantial but finite, while the magnetic 
correlations along the c axis are very short-range .13 It was 
difficult to investigate the in-plane correlation length at 
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FIG. 3. Scans of the (040) o /(400) o Bragg peaks at several 
temperatures through the transition from the LTO phase (a) , 
to the LTT phase (d), with coexistence of the two phases at 
intermediate temperatures, (b) and (c). Measured at H7 us- 
ing Et = 14.7 meV, horizontal collimations = 10'-10'-10'-10', 
monochromator = PG(002), analyzer = PG(004). 
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FIG. 4. Elastic scans in the transverse direction (see in- 
set) through the magnetic superlattice peak at (0.5,0.382) 
in Lai.48Ndo.4Sro.i2Cu04 for several different temperatures. 
The solid lines through the data points are fits to a Lorentzian 
plus constant background. Note the change in the inten- 
sity scale between (b) and (c). Measured at H4M using 
Er = 14.7 meV, horizontal collimations = 40'-80'-20'-80'. 
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FIG. 5. Results of the Lorentzian fits to the magnetic peak 



scans, as 



function of temperature, (a) Integrated intensity 
(in arbitrary units); arrow indicates the magnetic ordering 
temperature detected by /xSR (Ref. |22| . (b) Peak half- width, 
k; dashed line indicates resolution half-width; line through 
points is discussed in the text. In fitting the scan at T = 55 K 
(squares), the width was constrained to lie on the curve, and 
only the amplitude was varied. 



higher temperatures because the weakness of the inten- 
sity limited the degree to which it was practical to obtain 
sufficient Q resolution with tightened collimation. 

Making another attempt, we have now succeeded in fol- 
lowing the temperature dependence of the peak widths 
by performing transverse scans, where we can take ad- 
vantage of the intrinsically tight resolution at small Q. 
Examples of such scans are presented in Fig. |4|. The fit- 
ted peak shape corresponds to a Lorentzian plus a con- 
stant background. (No systematic temperature depen- 
dence was found for the background, so the background 
value averaged over essentially all temperatures was ap- 
plied and held fixed in the final fits.) An alternative peak 
shape has been suggested based on the recently proppsed 
analogy between charge stripes and liquid crystals.E3 If 
the pinned stripes behave like a smectic liquid crystal, 
then one might expect the spatial correlations to decay 
algebraically, rather than exponentially^ In that case, 
the peak shape in reciprocal space would be proportional 
to (Q — Qo) p - However, if the exponent p is close to 2, 



one would need a much better signal-to-noise ratio and 
excellent knowledge of the background in order to distin- 
guish the algebraic line shape from a Lorentzian. Further 
experimental advances will be required to test for such 
differences. 

Results for the integrated intensity and peak half- 
width at half-maximum, k, are shown as a function of 
temperature in Fig. |^. As indicated in (b), the resolu- 
tion width is much less than the measured width, thus 
justifying the absence of a resolution correction in the 
fitting. Considering first the intensity, it is interest- 
ing to compare with a recent muon-spin-rotation (/iSR) 
study of magnetic order in a similar crystal.E3 There, a 
loss of order was detected near 30 K, as indicated by 
the arrow in Fig. |5|(a). In contrast, the "elastic" signal 
seen by neutron scattering appears to survive up to ap- 
proximately 50 K.r-Such an apparent conflict has been 
found previouslyEjH in studies of magnetic correlations 
in La2- £C Sr, r CuC)4 with x ~ 0.05. It can be understood in 
terms of the distinct sensitivities of the two techniques to 
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fluctuations. In /iSR, the local hyperfine field is consid- 
ered static when it does not vary significantly on the time 
scale of the muon lifetime, ~ 2 fxs. On the other hand, the 
neutron measurement integrates over fluctuations within 
a Gaussian resolution window with a full width at half 
maximum of approximately 0.8 meV (0.2 THz). Thus, 
the observations indicate that between 30 K and 50 K 
there is no static order, but instead there is a charac- 
teristic fluctuation rate that is gradually increasing from 
10 6 to 10 11 Hz. 

The temperature dependence of k is qualitatively con- 
sistent with the above scenario. Below 30 K, k remains 
constant. Taking the inverse of n at low temperature 
gives a correlation length of 200 A. Above 30 K, k starts 
to increase. The curve through the data points corre- 
sponds to a fit with the form 



fitted value of B we find a rs-0.01. Going further, the 
theoretical expression for A iscJ 



K 



A, 



,-B/knT 



(1) 



which is essentially the formula used by Keimer et al. 
to describe the inverse correlation length in lightly doped 
La2Cu04. The fit yields the following parameter val- 
ues: k = 0.0050(2) A" 1 , A = 0.6(6) A" 1 , and B = 
17(4) meV, where the estimated uncertainty in the last 
digit is listed in parentheses. 

The second term on the right-hand side of Eq. (|l|) 
has the form predicted for the 2D quantum antiferro- 
magnetie JJeisenberg model in the renormalized classical 
regime.L3'E3 Fits to experimental measurements of k(T1 
in pure and lightly-doped La2Cu04 by Keimer et alSB 
yielded B pa J, where J is the exchange energy between 
nearest- neighbor Cu spins. It is not immediately clear 
that such a formula, based on a spin-only model, should 
be relevant to the present case, where there is a substan- 
tial cpacentration of holes in the planes. Castro Neto and 
Honec.3 have argued that within the stripe model, where 
the holes are restricted to moving along antiferromag- 
netic domain walls, the main effect of the holes would 
be to weaken the effective exchange between spins on ei- 
ther side of a domain wall. Assuming that this weakened 
coupling is the dominant factor for determining long- 
wavelength excitations along the modulation direction, 
they obtain an effective anisotropic Heisenberg model. 
Analysis of k(T) based on that model yields the sam 
exponential dependence on the inverse temperature^j'CJ 
as obtained in the isotropic case, but with modifications 
to the relationship between B and J. 

Let a represent the value of the effective exchange per- 
pendicular to the stripes relative to the exchange parallel 
to the stripes. Then, according to Ref. H 



with 



B = 2Trp s {a), 



p s (a) = y/aJS 2 



(2) 



(3) 



where S = | is the spin. If we assume J to have the 
same value as in La2Cu04 (pa 135 meVES), then from the 



A pa 8V2wa^/ae, 



(4) 



where a is the lattice parameter. Plugging in a pa 0.01 
gives A ps 0.5, in good agreement with the fitted-|¥alue. 

Within the anisotropic Heisenberg model,E-iE3 the 
spin- wave velocity perpendicular to the stripes (along the 
modulation direction) is 



2 C0 ' 



(5) 



where Co is the spin-wave velocity of the undoped anti- 
ferromagnet.—The fitted value of a, together with the 
experimental^ cq for La2Cu04, gives hc± pa 50 meV-A. 
We can use this result to check the assumption that the 
neutron measurement is integrating over low-energy fluc- 
tuations for T > 30 K. Neutrons measure the dynamic 
susceptibility, x"(Q, uj), multiplied by 1/(1 - e - fiw / feBT ). 
Integrating over Q, there is a peak in this quantity at 
r pa Tick. Using the numbers above, we find T pa 0.2 meV 
at 30 K, and T pa 0.8 meV at 50 K, where the signal is dis- 
appearing. The former value is smaller than the energy 
resolution half- width, while the latter is larger. Thus, our 
interpretation appears reasonable. On the other hand, 
we will show in the next section that Hc± pa 50 meV-A is 
much too small to be compatible with the observed low- 
energy inelastic scattering. It is also incompatible with 
the range_aL energies over which incommensurability is 
observedoEJ in La2_ a; Sr a ;Cu04 with x pa 0.15. 

One solution to this_dilcmma is suggested by the analy- 
sis of Tworzydlo et alE3 They have pointed out that when 
the magnetic domains become sufficiently narrow (small 
stripe spacing) one can view them as coupled spin lad- 
ders. For small a, one obtains qualitative differences de- 
pending on whether the spin ladders have an even or odd 
number of legs. (This distinction is absent in thq-Qciginal 
treatment of the anisotropip-Heisenberg model.E5Ea) In 
particular, Tworzydlo et alr3 have shown that for 2-leg 
ladders a quanturrwdisordered state occurs for a < 0.30. 
(For 4-leg laddersj£3 the transition is at a = 0.07.) In 
terms of the 2-leg ladder model, our observed tempera- 
ture dependence of k could be explained with a pa 0.35, 
which jszpuld also be consistent with relatively stiff spin 
waves .E3 | — . 

Tworzydlo et alEB have suggested that the 2-leg lad- 
der model could be appropriate at x = | if the hole 
stripes are centered on O rowsE3 in a manner such that 
they tie up the spins on the immediately adjacent rows 
of Cu. If this were the case, then tuning x away from | 
should introduce odd-leg ladders, which would-presum- 
ably enhance the tendflscy for magnetic order .E3 In con- 
trast, experimentstl'caEJ indicate that the highest mag- 
netic ordering temperature actually occurs at x ~ 0.12. 
Thus, while the spin-only models can explain several fea- 
tures of the measurements, questions still remain. It is 
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quite possible that the fluctuations of the charge stripes, 
which are ignored in spin-only models, play a significant 
role in determining the spin correlations. 

One other puzzle concerns the saturation of n below 
30 K. If the characteristic energy for critical spin fluctu- 
ations were 0.2 meV (or perhaps higher), then most of 
the spin fluctuations would be too fast to allow any muon 
precession to be observable. One possibility is that our 
analysis is too simple minded. Alternatively, it may be 
that the lack of temperature dependence in k at low T is 
associated with disorder in the charge-stripe array, rather 
than with the temperature dependence of the spin fluc- 
tuations. Analysis of such a disorder model is presented 
in section VII. 



V. INELASTIC MAGNETIC SCATTERING 

It is also interesting to characterize the spin dynam- 
ics directly through measurements of inelastic scatter- 
ing. Figure ^ shows results for x = 0.12 obtained by 
scanning through the incommensurate positions, along 
Q = (h, 0.5,0), with the energy transfer fixed at 2 meV 
[(a),(b),(c)] and at 3 meV [(d),(e),(f)]. A linear back- 
ground has been subtracted, corrections have been made 
for differences in resolution volume at the two energies, 
and the thermal factor has been divided out, so that the 
resulting signal is proportional to %"(Q,w). The scat- 
tering is peaked rather sharply about the incommensu- 
rate wave vectors defined by the elastic scans, with the 
Q- width limited by resolution at 50 K. The signal is rel- 
atively strong at 50 K, and gradually decreases as the 
temperature is raised. At 72 K, above the transition to 
the LTO phase, the signal is still finite at 3 meV. ThisJast 
result provides an important connection with studiesc3'Eil 
of spin fluctuations in La2- a; Sr. I ;Cu04, where the struc- 
ture is LTO. The LTT phase induced by Nd substitution 
appears to stabilize a static component of correlations 
that already exist dynamically in the LTO phase. 

Inelastic measurements have also been performed on 
an x = 0.15 sample. Results for the local susceptibil- 
ity, x"{oj) = J <1Q2dx"(Qj k-Oi a t two excitation ener- 
gies, 1.75 and 3.5 meV, are presented as a function of 
temperature in Fig. fj]. x"{uj) changes relatively little be- 
tween 10 and 40 K, where an elastic magnetic signal is 
observed. Above the neutron-determined ordering tem- 
perature (~ 45 K), the signal at 1.75 meV falls off rapidly, 
whereas the susceptibility at 3.5 meV appears to decrease 
more slowly. This behavior is qualitatively similar to 
the decay of the low^energy spin-fluctuation intensity in 
La2Cu04 above TjyE-i 

Finally, it is interesting to compare the Q-width of 
the excitations in the x = 0.12 crystal with that found 
in La2- x Sr x Cu04 with x = 0.15. Figure |^ shows scans 
for two such samples obtained under similar conditions: 
Tiui = 3 meV and T = 40 K. The curve through the data 
in (b) is a Gaussian plus a linear background, while the 



La i.4 8 Nd 0.4 Sr 0,2 CU °4 



Q = (/?, 0.5, 0) 



~i i i i r 
fvu> = 2 meV 



O 



O 



3 

a 




0.3 0.5 0.7 0.3 0.5 0.7 
h h 

FIG. 6. Constant-energy scans along Q = (h, 0.5, 0), 
through the magnetic wave vectors h — 0.5 ± e, for 
Lai.48Ndo.4Sro.i2Cu04. The Bose factor has been divided 
out, a fitted background (linear in h) has been subtracted, 
and a correction has been made for differences in resolution 
volume at the two energies. Panels on the left correspond to 
Tkij — 2 meV and temperatures of (a) 72 K, (b) 68 K, (c) 50 K; 
on the right, hu — 3 meV and temperatures are (d) 72 K, (e) 
68 K, (f) 50 K. Measured at H4M using a fixed incident en- 
ergy of 14.7 meV, horizontal collimations = 40'-120'-80'-80', 
and typical counting time of 20 min per point. Solid lines are 
gaussians fit with the constraint that the peaks be symmetric 
about h = 0.5. 
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FIG. 7. Q-integrated x"(Q> w ) for Lai.45Nd .4Sr .i5CuO4 
obtained from scans similar to those in the previous figure. 
The vertical line indicates the LTT-LTO phase boundary, and 
the dashed line is a guide to the eye. 



curve in (a) is a fit based on a disordered stripe model 
described in Ref. [l5|. The peak in (b), which is resolution 
limited, is clearly much narrower than that in (a). Such 
a difference appears to be consistent with the more dy- 
namic nature of the correlations in the x = 0.15 sample 
with no Nd. 

The inelastic Q width at 2 meV in La2_ x Sr x Cu04 
with = 0.12 has recently been reported by Yamada 
et aZ.EEl They find a value essentially the same as the 
resolution-limited width that we observe at 3 meV in the 
Nd-doped crystal. An elastic magnetic signal at the same 
incommensurate position haSj-also been discovered in the 
x = 0.12 crystal with no Nd.E3 Thus, there seems to be 
a reasonable correlation between decreasing inelastic Q- 
widths (perhaps limited by the correlation length) and 
the onset of static order. 

Before moving on, we should consider whether the ef- 
fective spin-wave velocity discussed in the last section 
is consistent with the observed peak width (half-width 
= 0.012 A" 1 ) in Fig. |(b). If the width were determined 
entirely by dispersion, rather than resolution, then the 
corresponding spin-wave velocity is ~ 250 meV-A. This 
is almost 3 times greater than the estimate we obtained 
for TiCeff from n(T). Part of the discrepancy might be as- 
sociated with anisotropy in the dispersion; however, at- 
tempts to test this experimentally have so far been lim- 
ited by anisotropy in the resolution function. Clearly, 
more work is required to obtain a consistent quantitative 
account of the quasi-elastic and inelastic scattering. 
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FIG. 8. Comparison of constant-energy 

scans at Tiuj = 3 meV through an incommensurate magnetic 
peak (see inset) for (a) Lai.85Sro.isCu04 (Ref. ^), and 
(b) Lai.48Ndo.4Sro.i2Cu04 (crystal U2). Both measurements 
were performed with a fixed incident energy of 14.7 meV and 
a sample temperature of 40 K. Scan (a) was done at H7 with 
horizontal collimations of 40'-40'-40'-80', while (b) was done 
at H4M with 40'-80'-80'-80'. Fitted curves are discussed in 
the text. 
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FIG. 9. Elastic scans along Q = (ft, 0, 2) at two different 
temperatures for Lai.45Ndo.4Sro.i5Cu04. Note that the in- 
tensity is plotted on a logarithmic scale. The strong peak at 
h m 2.975 corresponds to the LTO-allowed (032) superlattice 
peak. The broad feature centered near h = 3 at 80 K is dif- 
fuse scattering associated with the (302) reflection, which is 
forbidden in LTO but allowed in LTT and LTLO. The fitted 
curves are discussed in the text. 



VI. IMPERFECT STRUCTURAL ORDER IN THE 
LTO PHASE 

We now return our attention to the structural order of 
the x = 0.15 crystal in the LTO phase. In section III we 
noted the presence of a broad, unexpected peak between 
the (040) o and (400) o Bragg reflections [see Fig. |(a)]. 
The position of the broad peak suggests a tetragonal 
(or less-orthorhombic) component, and the width cor- 
responds to domains of width ~ 140 A. To test whether 
this component is consistent with LTT-like domains, we 
looked for related scattering near a characteristic super- 
lattice position. Figure || shows scans along Q = (h, 0, 2) 
(note the logarithmic intensity scale). The strong peak 
at h w 2.975 is the allowed (032) o of the LTO phase; the 
noninteger value of h is due to the slight orthorhombic 
splitting between a* and b*. In the 80 K scan, one can 
see a second peak, broad and weak, that is centered at 
h 1=3 3.00. Such a peak is not allowed in the LTO phase, 
but would be expected in the LTT. Its intensity has di- 
minished considerably in the scan at 230 K. 

To analyze the temperature dependence of the diffuse 
(302) o scattering, we restricted our choice of fitting func- 
tions to Gaussians. (Lorentzians were rejected because 
of the significant peak area associated with tails outside 
of the measurement region.) A single gaussian was used 
for the (032) o peak, but was found to give a poor fit to 
(302)o- A more satisfactory fit was obtained using the 
following combination of two Gaussians: 



1(h) = A 



-±(h-ll ) 2 /cr 2 



-±(h-h Q ) 2 /(2a) 2 



(6) 



Typical fits arc indicated by the curves through the data 



points in Fig. |. At 80 K, the peak center, ho, is 2.992(2), 
and the widths of the two Gaussian peaks correspond to 
domain sizes of 180 A and 90 A. The average of these 
lengths is similar to the value obtained from the broad 
(400) o peak. As the peak width does not vary signifi- 
cantly with temperature, the intensity reflects the volume 
fraction of LTT-like domains within the LTO matrix. 

The temperature dependence of the (302) o intensity is 
plotted in Fig. [l0](a). The intensity varies little at high 
temperature, but it begins to grow rapidly at lower tem- 
peratures. As indicated by the curve through the data 
points, the temperature dependence is described well by 
a function of the form a + b/T 3 . The significance of the 
temperature independent term is not clear. As one can 
see in the 230-K scan in Fig. |[ the diffuse peak becomes 
difficult to distinguish from possible (very weak) tails on 
the strong (032) o peak. Thus, it is possible that the 
scattering from LTT-like domains actually disappears at 
room temperature. 

It is interesting to compare the temperature depen- 
dence of the (302) o scattering with that of the in-plane 
resistivity. In the LTO phase, p a b is nearly linear in tem- 
perature, just as is observed for samples with no Nd; 
however, a slight deviation from linearity begins to be- 
come noticeable below 150 K. To examine this deviation, 
we first fit the data between 170 and 300 K to the form 

Pab = Pq + PtT, (7) 

finding po = —0.01 fiil-cm and px — 1-26 x 10 -3 /xf2- 
cm/K. Extrapolating the linear contribution and sub- 
tracting from the data gives Ap a t, which is plotted in 
Fig. |l(](b). As one can see, Ap a b grows at low tem- 
perature in a fashion similar to the (302) o intensity. In 
fact, the temperature dependence is well described by the 
same functional form, as indicated by the curve though 
the data points. 

Given the evident correlation between the (302) o in- 
tensity and Ap a h, we speculate that the rise in Ap a b 
is associated with charge-stripe pinning in the LTT-like 
domains which occupy a small volume fraction of the 
sample. Electron-diffraction studies of Lai.88Bao.i2Cu04 
and Lai.5Ndo.4Sro.iCu04 indicate that LTTJikc-domains 
first appear at the LTO twin boundaries B-HeI If LTT- 
like domains uniformly decorate all LTO twin bound- 
aries, then it may be difficult for an electron (or hole) 
to move through the sample without passing through a 
region of pinned or nearly-pinned charge stripes. The 
continuous development of the LTT-like volume fraction 
and its impact on the resistivity may explain, at least 
in part, why the resistivity does not clearly reflect the 
temperature dependence ftfj-'ihe charge order parameter 
determined by diffraction.l3'B3 



VII. SUMMARY AND DISCUSSION 

We have presented neutron scattering studies of vari- 
ous aspects of ordering in crystals of 
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FIG. 10. (a) Temperature dependence of the diffuse in- 
tensity centered at Q — (3,0,2), for Lai.4sNdo.4Sro.i5Cu04. 
Solid line is a fit to the form / = a+b/T 3 ; dashed line indicates 
the temperature-independent term, (b) In-plane resistivity 
remaining after subtraction of a dominant contribution linear 
in T, determined from a fit over the range 170 K < T < 300 K. 
The line through the data is a fit to the form Ap a & = a + b/T 3 . 
The inset shows the data and fit over a larger temperature 
range. 



Lai. 6 _ a; Ndo.4Sr 2; Cu04 with x = 0.12 and 0.15. The 
most significant results are the following: 1) At low tem- 
perature, where the Cu spins order in a nearly two- 
dimensional stripe structure, the magnetic correlation 
length remains finite. 2) The temperature dependence of 
the correlation length above 30 K is similar to that of the 
renormalized-classical regime of a 2D Heisenberg antifer- 
romagnet, and applying the formula for the spin-only sys- 
tem yields a measure of the effect that the charge stripes 
have on the effective spin-wave velocities. 3) Within the 
LTO phase, there is a correlation between the deviation 
of the resistivity from a linear temperature dependence 
and the volume fraction of LTT-likc domains. Below, 
each of these points will be discussed in turn. 

A. Finite correlation length 

What can we say about the nature of the magnetic dis- 
order in the ground state of the x = 0.12 sample? One 
possible model would be isolated, ordered domains of fi- 
nite size within a sea of disorder; however, such a picture 
is in conflict with the zero-field /iSR observation that es- 
sentially every muoH-implanted in such a sample sees a 
local hyperfine field£3 There appear to be no significant 
regions of the sample that do not show local magnetic 
order. An alternative model that allows local order ev- 
erywhere involves disorder in the spacing between charge 
stripes. One assumes that the charge stripes are sepa- 
rated by an integral number of lattice spacings, and that 
at least two distinct stripe periods exist. Such a model is 
consistent with the idea that commensurability with the 
lattice should play a role in the pinning of stripes in the 
LTT phase. Furthermore, as we will see, it is compatible 
with the incommensurate peak splitting,!!! e = 0.118. 

Formulas for calculating the scattered intensity from a 
random mixture of 2 structural— units were first worked 
out by Hendricks and Teller .l1!e3 (One of us used these 
formulas previously^ to model staging disorder associ- 
ated with interstitial oxygens in La2Ni04+,5.) As struc- 
tural units we take magnetic unit cells of period 8a and 
10a. (Each of these corresponds to two unit cells of the 
corresponding charge modulation.) Within each unit, we 
consider only Cu sites, and assume, for simplicity, a si- 
nusoidal variation of the spin density. The atomic dis- 
placements associated with the charge modulation are 
also taken to be sinusoidal. 

Within the model, the only free parameter is the ra- 
tio of the occurrence probabilities for the two structural 
units. We adjusted this value to make the peaks associ- 
ated with spin modulation appear at positions consistent 
with experiment. This criterion was satisfied with a prob- 
ability of 75% for the 8a unit and 25% for the 10a unit. 
The resulting scattered intensities due to the spin and 
lattice modulations are shown in Fig. |ll| as a function of 
momentum (in reciprocal lattice units) along the mod- 
ulation direction. To extract peak positions and widths 
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FIG. 11. Calculated scattered intensity as a function of h 
for the model of ID stripe-spacing disorder discussed in the 
text. Thick solid line: scattering due to spins; thick dashed 
line: signal due to lattice modulation. Thin lines represent 
fitted Lorentzians. 

TABLE I. Values for effective e (from peak position) 
and half-width at half maximum (HWHM) obtained from 
Lorentzian fits to the disordered-stripe model calculation 
shown in Fig. [ill 



Peak 


e e ff 


HWHM 




(rlu) 


(io- 3 A" 1 ) 


Spin 


0.118 


9.5 


Lattice 


0.120 


38.3 



for comparison with experiment, Lorentzians were fit to 
the calculated intensity. The results are listed in Table |. 

It is obvious from the figure and the fitted values in the 
table that the disorder assumed in the model yields finite 
peak widths that are larger for the lattice peaks than for 
the spin peaks. The width for magnetic peaks in this sim- 
ple model is about twice what we have found experimen- 
tally, and the. same discrepancy is found for the charge 
order peaks.o However, the ratio of calculated widths 
for lattice peaks relative to spin peaks (= 4) is exactly 
the sanaa as the ratio of the experimentally determined 
widths .E2I The difference in widths for lattice and spin 
peaks is easily understood: the scattering from the two 
distinct unit cells of the model is more nearly in phase at 
the magnetic wave vector than at the charge-order wave 
vector. While one could tune the choice of unit cells in 
the model to obtain better quantitative agreement with 
experiment, the salient point here is that the experimen- 
tally observed peak widths and positions are qualitatively 
consistent with disorder in the stripe spacing. 

What causes the stripe disorder? One possibility is 
the local potential variations due to the randomly posi- 
tioned divalent Sr ions. JDn the other hand, evidence has 
recently been prcsentcdL3 that all variations in cation ra- 
dius on the La site contribute to the suppression of T c , so 
perhaps the substituted Nd, which is significantly smaller 



than the La, makes a contribution to disorder compara- 
ble to that of the Sr. The relevance of disorder to pinning 
charge stripes and the development of glassy order is dis- 
cussed in Refs. [3f| |f| 

In applying the Hendricks- Teller model, we have con- 
sidered disorder in only one direction, perpendicular to 
the stripes. If the disorder is induced by randomly 
positioned Sr and Nd substituents, then there should 
also be defects, such as dislocations, along the stripes. 
Equal probabilities for defects in parallel and perpendic- 
ular directions would explain the absence of a detectable 
anisotropy in scattering peak width. 

B. Temperature-dependence of the correlation 
length 

We have seen that the magnetic peak width begins 
to grow for temperatures above 30 K, the point at which 
/iSR measurementsE3 indicate the disappearance of static 
magnetic order. Given the finite widths for both the 
magnetic and charge-order peaks at low temperature, it 
is reasonable to ask whether the decrease in the magnetic 
correlation might be due to increasing disorder (perhaps 
dynamic) in the positions of the charge stripes. Accord- 
ing to the Hendricks- Teller model applied above, such 
an increase in charge stripe disorder should be reflected 
in a relative growth in the charge-order peak width 
comparable to that of the magnetic peaks. However, 
the charge-order width observed in the X-ray diffraction 
measurements^ shows little variation up to 65 K. Hence, 
it appears most likely that the charge stripe configura- 
tion remains static throughout the region in which the 
magnetic width is observed to vary. 

It seems reasonable, then, to analyze the decrease of 
the magnetic correlation length in terms of spin fluctu- 
ations within well defined domains. In section IV we 
considered the implications of an effective anisotropic 
HeiseBbeEg model first proposed by Castro Neto and 
Hone.E5c3 The temperature dependence of k is well de- 
scribed by the renormalized-classical behavior predicted 
by the analysis of the model; however, the effective spin- 
wave velocity perpendicular to the stripes implied by 
the fit to the model appears to be inconsistent with the 
measurements of the spin fluctuations at tuv = 3 meV. 
The model can be improved by considering the enhanced 
quantum effects that occur when the magnetic domains 
correspond to even- legged spin ladders. The latter model 
appears to run into trouble when one tries to consider the 
doping dependence of the maanetic correlations, since 
magnetic order is most robustlilleSli-l at x » i. 

It is quite possible that the interaction between the 
holes in the domain walls and the spins in the domains is 
not adequately characterized in terms of an effective ex- 
change coupling between antiferromagnetic domains. An 
alternative model of the interaction has been evaluated 
by Emery, Kivelson, and ZacharB They find that the in- 
teraction takes place by the hopping of pairs of holes from 
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a charge stripe into the magnetic regions, and that the 
pair hopping tends to enhance singlet correlations. This 
picture would certainly be consistent with a weak effec- 
tive coupling between neighboring domains; however, a 
specific prediction for the low-energy spin fluctuations 
relevant to the present case has not been made. 

C. Correlation between resistivity and structure 

For crystals of Lai.g_ a; Ndq^Sr x Cu04 with x near g, 
the in-plane resistivity showsLJ an abrupt upward jump 
on transforming to the LTT structure near 70 K. When 
Nd is replaced by Eu, the LTT transition shifts up to ap* 
proximately 130 K, and no jump in resistivity is found.tij 
Instead, on cooling there is a continuous deviation from 
a linear temperature dependence, with a significant up- 
turn below 50 to 60 K. In light of that, the correlation 
between resistivity and LTT-like volume fraction above 
80 K in our Nd-doped x = 0.15 sample should perhaps 
come as no surprise. 

It seems reasonable to assume that the rise in resis- 
tivity on cooling is associated with pinning of the charge 
stripes. If the charge order were like that in a conven- 
tional CDW system, in which a gap opens over a signifi- 
cant portion of the Fermi surface, then one might expect 
to see the charge order parameter directly reflected in 
the resistivity. So far, infrared reflectivity studies have 
provided no evidence for a charge gap JIB and such a gap 
would appear to be incompatible with jtte-SHneiconduc- 
tivity observed in the x — 0.15 sample JilE2lOE3 On the 
other hand, the gradual increase in resistivity is con- 
sistent with the glass-like stripe disorder that was an- 
alyzed above. Also, a recent study of the CDW system 
TaS2 has shown that a very small amount of disorder in 
that system can suppress the transition to commensuxate 
CDW order and the concomitant jump in resistivity^ It 
has been argued recently that transverse fluctuations of 
charge stripes are important for suppressing CDW or- 
der along the stripes.Ea Perhaps static disorder caused 
by randomly positioned dopants can also frustrate true 
CDW order. 
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